Uncertainties in modeling low-energy neutrino induced reactions on iron group nuclei 
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Charged-current neutrino-nucleus cross sections for ^''-^^Fe and ^*'^"Ni are calculated and com- 
pared using frameworks based on relativistic and Skyrme energy density functionals, and the shell 
model. The current theoretical uncertainties in modeling neutrino-nucleus cross sections are as- 
sessed in relation to the predicted Gamow- Teller transition strength and available data, multipole 
decomposition of the cross sections, and cross sections averaged over the Michel flux and Fermi-Dirac 
distribution. Employing different microscopic approaches and models, the DAR neutrino-^^Fe cross 
section and its theoretical uncertainty are estimated: < a >th = (258±57) xlO~^^cm^, in very good 
agreement with the experimental value: < a >ea;p = (256±108±43) x lO^^^cm^. 

PACS numbers: 21.30.Fe, 21.60.Jz, 23.40.Bw, 25.30.-c 



Weak neutrino-induced processes in nuclei provide in- 
formation of relevance for modeling the response in neu- 
trino detectors, understanding the fundamental proper- 
ties of the weak interaction, and the role of neutrinos in 
stellar environment. Data on neutrino- nucleus cross sec- 
tions are presently available only for -'^^C and ^®Fe target 
nuclei, obtained by the LSND [i| and KARMEN [1 1| 
collaborations, and at LAMPF ]j|. At present only the- 
oretical approaches can provide cross sections for a large 
number of target nuclei that are involved in various ap- 
plications of neutrino physics and astrophysics. It is, 
therefore, crucial to quantitatively assess the theoretical 
uncertainties in modeling neutrino induced processes, in- 
cluding the detailed structure of principal transitions in- 
volved, and the total cross sections averaged over selected 
neutrino fluxes. The evaluation of current theoretical un- 
certainties in modeling neutrino induced processes is also 
important in view of future experimental programs, e.g., 
spallation neutron source (SNS) at ORNL Q, Large Vol- 
ume Detector in Gran Sasso (LVD) Q, and the beta- 
beams for the production of neutrinos by using the /3- 
decay of boosted radioactive ions [3, Q . 

Over the years a variety of microscopic models have 
been developed and employed in the calculation of 
neutrino-nucleus cross sections at low energies. These 
include the shell model [g-[l3| , the random phase approxi- 
mation (RPA) [il,[il,[il, continuum RPA (CRPA) Kdh 
[Tsj . hybrid models of CRPA and shell model fl4 . fl9| . 
Fermi gas model [23|, quasiparticle RPA [2l|, [22], pro- 
jected QRPA [2^, and relativistic quasiparticle RPA 
(RQRPA) [23 ■ For the purpose of the present analysis 
of theoretical uncertainties in modeling neutrino-nucleus 
cross sections, we chose the iron group nuclei for which 
the framework based on the energy density functionals 
and the shell model represent feasible approaches and, 
in addition, data from muon decay at rest (DAR) are 
available [3, [3| ■ The framework based on energy density 
functionals employs the Relativistic Hartree-Bogoliubov 
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FIG. 1. (Color onhne) GT transition strength in ^''Fe, cal- 
culated using the RQRPA (DD-ME2) and the Shell model 
(GXPFU). 



model (RHB) to determine the nuclear ground state, 
and the RQRPA to calculate all relevant transitions in- 
duced by the incoming neutrinos 2J] . Model calculations 
are performed using effective interactions with density- 
dependent meson- nucleon couplings, in this case the DD- 
ME2 interaction ^2^ , whereas pairing correlations are de- 
scribed by the finite range Gogny force [2g| . The nuclear 
shell model employed in the present study is based on the 
GXPFU effective interaction ^ for A'' = 1+ channel, 
supplemented by the RPA based on a Skyrme functional 
(SGII) for other multipoles [l3|. A detailed analysis of 
(anti)neutrino - ^^Fe cross sections based on the QRPA 
with Skyrme functionals is given in Ref. [21[ . Shell-model 
calculations are carried out with the code MSHELL um . 



Because of its importance for neutrino-nucleus cross 
sections at low energies, we start by analyzing the 
Gamow- Teller(GT) transition strength in the iron group 
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TABLE I. GT transition strengths calculated using the 
Skyrme RPA (SGII, SLy5), RQRPA (DD-ME2), and shell 
model (GXPFIJ) jial f all using g^ = 0.743a), compared to 
experimental values [29l |. 
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FIG. 2. (Color online) GT transition strengths calculated 
with the Skyrme RPA (SGII, SLy5) and RQRPA (DD-ME2), 
in comparison to the shell model (GXPFIJ) [Q and experi- 
mental values |23|. (Q)RPA calculations include the quench- 
ing factor 0.8 in the axial- vector coupling constant qa- 
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FIG. 3. (Color online) Same as described in the caption to 
Fig[2] but for the GT"'' transitions. Experimental values are 
from Refs. ilH. 



nuclei. Fig. [T] displays the GT strength distributions for 
^^Fe, obtained using the sheh model (GXPFIJ) and the 
RQRPA (DD-ME2). In both cases the calculated tran- 
sition strength is folded by a Lorentzian with the width 
r=0.5 MeV. The RQRPA includes only 2qp configura- 
tions and, therefore, cannot provide the detailed struc- 
ture of excitation spectra obtained by the shell model. 
Nevertheless, one can observe that the calculated transi- 
tion strength distributions are in reasonable agreement. 
In Figs. [2] and |3] the GT='= transition strengths con- 
tributing to the neutrino-induced processes are compared 
for a set of iron group nuclei: 54,56pg g^^j 58,60]\jj rpj^g 

following models and respective parameterizations have 
been used: i) RPA based on Skyrme functionals (SGII, 
SLy5), ii) RQRPA (DD-ME2) ^25], and iii) shell model 
(GXPFIJ) ^]. Results of model calculations are 



compared with the data for GT" [291 and GT+ [30|-l33 
transition strengths. The (Q)RPA calculations include 
the quenching of the free-nucleon axial-vector coupling 
constant gA = 1.262 -^ gA = 1, corresponding to the 
quenching factor 0.8 in the GT transition operator. A 
quenching factor is also used in the shell model. How- 
ever, its value 0.74 is adapted to the effective interaction 
and model space under consideration 13]. One notices 
in Figs. [2] and [3] that the shell model reproduces the ex- 
perimental values of GT~ transition strength with high 
accuracy (except for ^°Ni), and also the GT"*" strength 
is reasonably reproduced. (Q)RPA based approaches, 
however, even by quenching the value gA ^ ^, overes- 
timate both the GT~ and GT"*" transition strengths. For 
^"''^^Fe and ^^Ni the relativistic QRPA resuhs for the 
B(GT^) are within experimental uncertainties. The fact 
that QRPA calculations systematically overestimate the 
measured GT* transition strength, even though different 
effective interactions are used, indicates that a somewhat 
stronger quenching of the axial- vector coupling constant 
might be necessary to reproduce the data. Actually, if 
the same quenching factor 0.74 used by the shell model 
is also employed in (Q)RPA calculations, a very good 
agreement with the shell-model results is obtained. The 
(Q)RPA and shell model B(GT~) values, all obtained us- 
ing the same quenching factor 0.74, are shown in Table 
I in comparison to the data [2^ . In this case for all four 
nuclei: 54,56pg ^^^^ 58,60]Xfj^ ^i^g (Q)RPA B(GT-) values 
are found in good agreement with the shell-model results, 
particularly for the relativistic QRPA (DD-ME2). We 
have verified that the same result is also obtained for the 
B(GT+) channel. This similarity is not obvious as the 
two theoretical frameworks have different foundations, 
use different effective interactions and model spaces, that 
is, the underlying structure cannot be directly compared, 
except for the resulting observables that one confronts to 
data. This result should also be considered with caution 
because of well known problem of missing GT strength, 
due either to excitations that involve complex configu- 
rations at higher excitation energies, or excitations that 
include non-nucleonic degrees of freedom. 

As already emphasized in previous studies of neutrino- 
nucleus reactions [2l|, [2J], not only GT^ transitions but 



g 1 E^=40 MeV 
J 6- 

2 4- 

c 

2 - 

— 

^ 8 - E^=6() MeV 

"l 6- 

'2 4- 
e 

2 - 

— 

8 ^ E =80 MeV 



pe ■ DD-ME2/g^ =C''-'?A 

Jgf DD-ME2 / g^* = 0.8g^ 
!■ GXPFIJ+SGII / g^'= 0.74g, 



fc^li li Ml 



I all isli 



FIG. 4. (Color online) Contributions of the multipole tran- 
sitions A'^ = — 4 to the inclusive cross section for the 
^''Fe(i^,,e~)^''Co reaction, at E^^ = 40, 60 and 80 MeV. Cal- 
culations include RHB-fRQRPA (DD-ME2), and the sheU 
model (GXPFIJ) (for the 1+ transition) plus the RPA (SGII) 
for higher multipoles. The quenching factors in qa are de- 
noted in the figure. 
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FIG. 5. (Color online) Inclusive neutrino-nucleus cross sec- 
tions in the 1+ channel for the ^^.sepg ^^^ ss.eojs^j target nu- 
clei, at the incoming neutrino energies E^^ = 40, 60 and 80 
MeV. The resuhs are obtained using the RHB-I-RQRPA (DD- 
ME2) framework, and the shell model (GXPFIJ). The axial- 
vector coupling g\ includes a quenching factor as denoted in 
the figure. 



also excitations of higher multipoles must be included, 
depending on the energy range under consideration. In 
this work we analyze the reaction ^^Fe(i/e, e~)^^Co in 
more detail, using two representative theoretical ap- 
proaches: the shell model (GXPFIJ) supplemented by 
the RPA based on Skyrme functionals (SGII), and 
the fully consistent relativistic framework RHB+QRPA 
(DD-ME2). The goal is to provide an estimate of theo- 
retical uncertainties of contributions from different mul- 
tipole transitions to the neutrino cross section. In Fig. |3] 
we plot the contributions of the multipole transitions 
A'^ = 0^ — 4^ to the inclusive cross section for the 
^'^Fe(i^e,e")'^^Co reaction, at E^^ = 40, 60 and 80 MeV. 
The RHB+RQRPA (DD-ME2) calculations include the 
standard (0.8), and enhanced (0.74) quenching factors 
of the axial-vector coupling constant gA in all multipole 
operators. In the non-relativistic framework the shell 
model (GXPFIJ) is used for 1+ transitions, and RPA 
(SGII) is used for higher multipoles. The quenching fac- 
tor for the axial- vector coupling qa is 0.74. In the shell 
model calculation of 1+ transitions the effect of finite 
momentum transfer (q) is taken into account by evaluat- 
ing the matrix elements < f\\jo{<if)\Y^ x t?]"^i_||i > and 
< f\\J2{qf')\Y'^ X (T]^i_||z > at each g, instead of the ap- 
proximate treatment of Ref. [13| where the GT matrix 
element < f\\ja[qR)at^\\i > was evaluated multiplied by 
jo{qR) (R is the nuclear radius). 

At relatively low neutrino energies {Ei, < 40 MeV) 
the dominant contribution to the calculated cross sec- 
tions originates from GT transitions (A'^=l"'"). With in- 
creasing E^, however, contributions from other multipole 
transitions become important. In particular, at Ei,=80 
MeV the dominant transition is the spin-dipole A'^=l~, 



but also other components, e.g. A'^ = l+, 2^, 2+, 3+, play 
an important role. The cross sections plotted in Fig. |3] 
show that the two models predict a very similar struc- 
ture and distribution of the relative contributions from 
various multipoles. These results are also in agreement 
with those discussed in Ref. |2l| . 

In Fig. [S] we compare the neutrino-capture cross sec- 
tions for the A'^^1"'" channel on the set of target nuclei: 
54,56pg g^j-^^ 58,60]\jj^ g^j^^j f^j. ^]-^g incoming neutrino en- 
ergies E„^ = 40, 60 and 80 MeV. The results are ob- 
tained using the RHB+RQRPA (DD-ME2) and the shell 
model (GXPFIJ). The axial- vector coupling 5^ includes 
a quenching factor as denoted in the figure. The cross sec- 
tions increase in heavier isotopes because electron neutri- 
nos are captured on neutrons. The two models, although 
based on different microscopic pictures, predict rather 
similar cross sections. One notes that at higher neutrino 
energies the shell model cross sections are slightly larger 
than those calculated with the RQRPA. Related to the 
previous discussion on the overall GT strength and the 
quenching of gA , from the cross sections shown in Fig. [5] 
it appears that the (Q)RPA does not require a stronger 
quenching than the usual 17^ = 1 to be in agreement with 
the shell model. The reason is that the calculated cross 
sections are not only determined by the overall GT tran- 
sition strength, but also by the transition energies that 
govern the energies of outgoing electrons. 

We have also analyzed cross sections averaged over the 
neutrino flux described by the Fermi-Dirac spectrum [2J] . 
Figure [6] displays the flux-averaged cross sections for the 
reaction ^^Fe{i'e, e~)^^Co, evaluated at different temper- 
atures in the interval T = 2 — 10 MeV, and for the 
chemical potential a=0. The RHB-I-RQRPA (DD-ME2) 
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FIG. 6. (Color online) Neutrino - ^^Fe cross sections aver- 
aged over the Fermi-Dirac distribution. The RHB+RQRPA 
(DD-ME2) calculations including the standard (0.8), and en- 
hanced (0.74) quenching factors of the a:xial-vector coupling 
gA , are compared to shell model -I- RPA results [13] , and those 
obtained using a hybrid model of Ref. [33 ]. 



calculations including the standard (0.8), and enhanced 
(0.74) quenching factors in gA, are compared to the shell 
model -I- RPA (GXPFIJ + SGII) results fll, and those 
obtained using a hybrid model of Ref. i33j. The latter 
model predicts somewhat larger cross sections, whereas 
a very good agreement is found between the results of 
shell model -I- RPA and RQRPA. 

The theoretical cross sections for the reaction 
^^Fe(i'e, e~)^^Co can also be analyzed in compari- 
son to data from the KARMEN collaboration. The 
calculated cross sections are averaged over the neu- 
trino flux described by the Michel spectrum ob- 
tained from muon decay-at-rest (DAR): f{E,y^) — 
(96i?^__/r7i^) {mfj. — '^E^^). Taking into account the re- 
sults obtained in this work with the RHB-f RQRPA (DD- 
ME2) - 263 X 10-42cni2 , shell model (GXPFIJ) (for 
1+ transitions) plus the RPA (SGII) for other multi- 
poles - 259 X 10~^^cm^ , as well as results from pre- 
vious studies that used the RPA with a Landau-Migdal 
force - 240 x lO-^^^m^ [H, the QRPA(SIII) - 352 x 



10- __ 

ism - 173.5 x lO^'^^cm^ [22]], the DAR neutrino-nucleus 
cross section and its theoretical uncertainty are estimated 
to be: < a >t,i=(258±57) xlO-^2f,jj^2^ r^j^jg ^^^^^ jg -^^ 

very good agreement with the data from the KARMEN 
collaboration: < a >ea;p=(256±108±43) x lO^'^^cm^. 
We note that the various models used to obtain the the- 
oretical estimate employ different effective interactions, 
and also comprise a wide range of values for the axial- 
vector coupling, from those without quenching [2l|,[22| to 
models that use a quenching factor of 0.7 [3J|. All the- 
ory frameworks, except the QRPA based on G-matrix 
formalism, favor the quenching of the axial-vector cou- 
pling constant gA, in accordance to constraints given by 
the experimental data on Gamow- Teller transitions. The 
implementation of the quenching of gA in QRPA based 
on Skyrme functionals 21] would lower the calculated 
neutrino-nucleus cross sections and it would further re- 
duce the overall theoretical uncertainty in < cr >th- 

In conclusion, the charged current neutrino- nucleus 
cross sections for ^4,56pg ^^i^ ss.eojsjj have been ana- 
lyzed by employing models based on the relativistic and 
Skyrme energy density functionals, and the shell model. 
The theoretical uncertainties in modeling neutrino in- 
duced processes have been examined by considering the 
Gamow- Teller transition strength and available data, the 
multipole decomposition of the calculated cross sections, 
and cross sections averaged over the Michel flux and 
Fermi-Dirac distribution. It has been shown that vari- 
ous models predict very similar multipole distributions of 
neutrino-nucleus cross sections. The corresponding cross 
sections averaged over the DAR neutrino spectra show 
that the current theoretical uncertainty, despite a variety 
of models and effective interactions that have been used 
in many studies, is actually smaller than the experimen- 
tal one, and could be even further reduced by constrain- 
ing the quenching of the axial-vector coupling constant 
gA to data for the Gamow- Teller transition strength. 
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